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ABSTRACT 

We present full sky microwave maps in five frequency bands (23 to 94 GHz) from the WMAP first 
year sky survey. Calibration errors are <0.5% and the low systematic error level is well specified. The 
cosmic microwave background (CMB) is separated from the foregrounds using multifrequency data. The 
sky maps are consistent with the 7° full- width at half-maximum (FWHM) Cosmic Background Explorer 
( COBE ) maps. We report more precise, but consistent, dipole and quadrupole values. 

The CMB anisotropy obeys Gaussian statistics with -58 < fsi < 134 (95% CL). The 2 < / < 900 
anisotropy power spectrum is cosmic variance limited for / < 354 with a signal-to-noise ratio >1 per 
mode to / = 658. The temperature-polarization cross-power spectrum reveals both acoustic features and 
a large angle correlation from reionization. The optical depth of reionization is r = 0.17 ± 0.04, which 
implies a reionization epoch of t r = 1 801|q° Myr (95% CL) after the Big Bang at a redshift of z r = 20 +J 0 
(95% CL) for a range of ionization scenarios. This early reionization is incompatible with the presence of 
a significant warm dark matter density. 

A best-fit cosmological model to the CMB and other measures of large scale structure works re- 
markably well with only a few parameters. The age of the best-fit universe is to = 13.7 zb 0.2 Gyr old. 
Decoupling was t dec = 379+* kyr after the Big Bang at a redshift of z d ec = 1089 ± 1. The thickness of the 
decoupling surface was A idee = 195 ± 2. The matter density of the universe is = 0 . 1351 q ; oo 9 ’ 
baryon density is = 0.0224 ± 0.0009, and the total mass-energy of the universe is f l toJ = 1 .02 ± 0.02. 

There is progressively less fluctuation power on smaller scales, from WMAP to fine scale CMB mea- 
surements to galaxies and finally to the Ly-a forest. This is accounted for with a running spectral index, 
significant at the ~ 2o level. The spectral index of scalar fluctuations is fit as n s = 0.93 ± 0.03 at wavenum- 
ber &o = 0.05 MpcT 1 (J € ff « 700), with a slope of dn s /d\nk = -0.03 l^jg in the best-fit model. This flat 
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universe model is composed of 4.4% baryons, 22% dark matter and 73% dark energy. The dark energy 
equation of state is limited to w < -0.78 (95% CL). 

Inflation theory is supported with n s « 1, Q tot « 1, Gaussian random phases of the CMB anisotropy, 
and superhorizon fluctuations implied by the TE anticorrelations at decoupling. An admixture of isocur- 
vature modes does not improve the fit. The tensor-to-scalar ratio is r(Jco - 0.002 Mpc" 1 ) < 0.71 (95% CL). 
The lack of CMB fluctuation power on the largest angular scales reported by COBE and confirmed by 
WMAP is intriguing. WMAP continues to operate, so results will improve. 

Subject headings: cosmic microwave background, cosmology: observations, early universe, dark matter, 
space vehicles, space vehicles: instruments, instrumentation: detectors, telescopes 


I. INTRODUCTION 

The cosmic microwave background (CMB) radiation was first detected by Penzias & Wilson (1965). After its 
discovery, a small number of experimentalists worked for years to better characterize the CMB as they searched for 
temperature fluctuations. A leader of this effort, and of the WMAP effort, was our recently deceased colleague, 
Professor David T. Wilkinson of Princeton University. He was also a leading member of the Cosmic Background 
Explorer (COBE) mission team, which accurately characterized the spectrum of the CMB (Mather et al. 1990, 1999) 
and first discovered anisotropy (Smoot et al. 1992; Bennett et al. 1992; Kogut et al. 1992; Wright et al. 1992). The 
MAP was recently renamed WMAP in his honor. 

The general recognition that the CMB is a primary tool for determining the global properties, content, and history 
of the universe has led to the tremendous interest and growth of the field. In addition to the characterization of the 
large-scale anisotropy results from COBE (Bennett et al. 1996; Hinshaw et al. 1996a, b; Kogut et al. 1996b, c, a; Gorski 
et al. 1996; Wright et al. 1996a), a host of experiments have measured the finer scale anisotropy (Benoit et al. 2003; 
Grainge et al. 2003; Pearson et al. 2002; Ruhl et al. 2003; Kuo et al. 2002; Dawson et al. 2001; Halverson et al. 2002; 
Hanany et al. 2000; Leitch et al. 2000; Wilson et al. 2000; Padin et al. 2001 ; Romeo et al. 2001 ; Harrison et al. 2000; 
Peterson et al. 2000; Baker et al. 1999; Coble et al. 1999; Dicker et al. 1999; Miller et al. 1999; de Oliverira-Costa 
et al. 1998; Cheng et al. 1997; Hancock et al. 1997; Netterfield et al. 1997; Piccirillo et al. 1997; Tucker et al. 1997; 
Gundersen et al. 1995; de Bemardis et al. 1994; Ganga et al. 1993; Myers et al. 1993; Tucker et al. 1993). As a 
result of these tremendous efforts, the first acoustic peak of the anisotropy power spectrum has been unambiguously 
detected (Knox & Page 2000; Mauskopf et al. 2000; Miller et al. 1999) and CMB observations have placed important 
constraints on cosmological models. Recently, Kovac et al. (2002) reported the first detection of CMB polarization 
arising from the anisotropic scattering of CMB photons at decoupling, ushering in a new era of CMB polarization 
measurements. 

The WMAP mission was designed to advance observational cosmology by making full sky CMB maps with 
accuracy, precision, and reliability, as described by Bennett et al. (2003a). The instrument observes the temperature 
difference between two directions (as did COBE) using two nearly identical sets of optics (Page et al. 2003c, a). These 
optics focus radiation into horns (Barnes et al. 2002) that feed differential microwave radiometers (Jarosik et al. 2003a). 
We produce full sky maps in five frequency bands from the radiometer data of temperature differences measured over 
the full sky. A CMB map is the most compact representation of CMB anisotropy without loss of information. 

In this paper we present the maps, their properties, and a synopsis of the basic results of the first-year of obser- 
vations. In §2 we give a brief overview of the WMAP mission. In §3 we summarize the data analysis, calibration, and 
systematic errors of the experiment, which are discussed in much greater detail in the companion papers by Hinshaw 
et al. (2003a), Page et al. (2003a), Jarosik et al. (2003b), and Barnes et al. (2003). In §4 we present the maps and their 
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sampling properties, and we compare the WMAP and COBE maps. In §5 we summarize the foreground analyses of 
Bennett et al. (2003a). In §6 we establish the Gaussian nature of the WMAP anisotropy, determined in the companion 
paper of Komatsu et al. (2003). In §7 we present the dipole and quadrupole moments and summarize analyses of 
the angular power spectrum (Hinshaw et al. 2003b; Verde et al. 2003). In §8 we highlight the WMAP polarization 
results, including a detection of the reionization of the universe (Kogut et al. 2003). In §9 we summarize some of the 
cosmological implications of the WMAP results (Page et al. 2003c; Spergel et al. 2003; Peiris et al. 2003). Finally, in 
§10 we discuss the availability of the WMAP data products. 


2. OBSERVATIONS 

The 840 kg WMAP Observatory was launched aboard a Delta II 7425-10 rocket (Delta launch number 286) on 30 
June 2001 at 3:46:46.183 EDT from Cape Canaveral. WMAP executed three phasing loops in the Earth-Moon system 
before a lunar-gravity-assist swing-by, a month after launch, catapulted WMAP to an orbit about the second Lagrange 
point of the Sun-Earth system, L2. Station-keeping is performed approximately four times per year to maintain the 
Observatory in a Lissajous orbit about the L2 point with the Earth- WMAP vector within about ~ 1° - 10° of the Sun- 
Earth vector. The phasing loop maneuvers and station-keeping are executed using the WMAP propulsion system of 
blow-down hydrazine and eight thrusters. 

The central design philosophy of the WMAP mission was to minimize sources of systematic measurement errors 
(Bennett et al. 2003a). The COBE mission proved the effectiveness of a differential design in minimizing systematic 
errors. Therefore, the WMAP instrument was designed with a back-to-back optical system with 1 .4 m x 1.6 m primary 
reflectors to provide for differential measurements of the sky. The primary and secondary reflectors direct radiation 
into two focal planes, with ten feed horns in each, as described by Page et al. (2003c). 

The beams have a gain pattern, G, which is neither symmetric nor Gaussian. We define the beam solid angle as 
f G(Q)/G max dQ. The beam size can be expressed as the square root of the beam solid angles, giving 0.22°, 0.35°, 
0.51°, 0.66°, and 0.88° for W-band though K-band, respectively. Alternately, the beams can be expressed in terms of 
a full- width at half-maximum (FWHM) for each band, given in Table 1. Detailed analyses of the WMAP beams are 
discussed by Page et al, (2003a, c). 

The feed horns are attached to orthomode transducers (OMTs) that split the polarization of the incoming signal 


Table 1 . Approximate Observational Properties by Band 


Item 

K-Band 

Ka-Band 

Q-Band 

V-Band 

W-Band 

Wavelength, A (mm) 

13 

9.1 

7.3 

4.9 

3.2 

Frequency, v (GHz) 

22.8 

33.0 

40.7 

60.8 

93.5 

Ant./therm. conversion factor, AT/AT a 

1.014 

1.029 

1.044 

1.100 

1.251 

Noise, (T 0 (mK) a = cr 0 AQj/ 2 

1.424 

1.449 

2.211 

3.112 

6.498 

Beam width 0(°FWHM) 

0.82 

0.62 

0.49 

0.33 

0.21 

No. of Differencing Assemblies 

1 

1 

2 

2 

4 

No. of Radiometers 

2 

2 

4 

4 

8 

No. of Channels 

4 

4 

8 

8 

16 
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into a differential correlation radiometer system with High Electron Mobility Transistor (HEMT) amplifiers. There 
are ten “differencing assemblies” each consisting of two “radiometers” with two “channels” each (Jarosik et al. 2003a; 
Bennett et al. 2003a). There are four W-band (~ 94 GHz), two V-band (~61 GHz), two Q-band (~ 41 GHz), one Ka- 
band (~ 33 GHz), and one K-band (~ 23 GHz) differencing assemblies. We usually refer to these bands by the generic 
designations K, Ka, Q, V, and W because there are multiple radiometers in each band, whose precise frequencies are 
not identical. Also, the effective frequency of a radiometer depends on the spectrum of the emission it detects. Precise 
frequencies for the radiometers for a CMB anisotropy spectrum are given by Jarosik et al. (2003a). Polynomials are 
given to determine the effective frequency of the radiometers depending on the emission frequency spectrum. See 
Table 1 for a summary of radiometer properties. 

Undesirable 1 // noise is minimized by the design of the WMAP radiometers (Jarosik et al. 2003a). All radiome- 
ters have 1// knees below 50 mHz; 18 of 20 are below 10 mHz; and 10 of the 20 are below 1 mHz Jarosik et al. 
(2003b). (The 1 // knee is defined as the frequency where the noise power spectral density is y/2 times higher than 
its high frequency value.) Jarosik et al. (2003a) demonstrate that all radiometer outputs have Gaussian noise, which 
“integrates down” with time as expected. 

The radiometers are passively cooled to ~ 90 K with no mechanical refrigerators. In addition, no actively cycling 
heaters were permitted anywhere on the WMAP spacecraft. These design features helped to ensure a mechanically, 
thermally, and electronically quiet platform that minimizes the driving forces of systematic measurement errors. 

In addition to the differential design, the COBE mission also demonstrated the importance of scanning large 
areas of the sky in a short period of time with a complex scan pattern. WMAP follows the COBE example with a 
three-axis (three reaction wheel) control system that maintains the Observatory in a nearly constant survey mode of 
operations. (The Observatory is in constant survey mode, except for only ~ 1 hr for each of ~ 4 station-keeping 
maneuvers per year.) In survey mode, the optical boresight sweeps out a complex pattern on the sky (Bennett et al. 
2003a). Approximately 30% of the sky is observed each hour. The Observatory spins at 0.464 rpm (~7.57 mHz) and 
precesses at 1 rev hr' 1 (~ 0.3 mHz). 

Six months are required for L 2 to orbit half way around the Sun, allowing for full sky coverage. The observations 
presented in this and companion papers include a full orbit about the Sun, thus containing two sets of full sky obser- 
vations. By 10 August 2001, WMAP was sufficiently stable in its L 2 orbit for CMB data-taking to commence. One 
year of observations, completed on 9 August 2002, were analyzed. Data taken beyond this date will be the subject of 
future analyses. 


3. DATA PROCESSING, CALIBRATION, & SYSTEMATIC ERRORS 

Time-ordered telemetry data from the Observatory are down-linked via NASA’s Deep Space Network (DSN) to 
the WMAP Science and Mission Operations Center (SMOC) at the Goddard Space Flight Center. The data are then 
transferred to the WMAP Science Team for analysis. All of the instrument data are down-linked to the ground without 
any on-board flight data processing, thus allowing full insight into potential systematic effects. 

Only a fraction of a percent of data was lost in the flow from the Observatory to the SMOC. About 1% of the 
received data were not used due to systematic error concerns (e.g., data taken during near station-keeping maneuvers). 
Of the ~ 99% good data, the processing pipeline flagged observations where bright planets were in the beams so that 
these data would not used be used in making maps. The statistics on lost, bad, and flagged data are given in Table 2. 

An overview of the data flow is shown in Figure 1 . The heart of the data analysis efforts center on studies 
of systematic measurement errors (Hinshaw et al. 2003b). Components of spurious signals at the spin period are 
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Fig. 1 . — An overview of the WMAP data flow. The references are: ( 1 ) Barnes et al. (2003) (2) Bennett et al. (2003b) 
(3) Bennett et al. (2003c) (4) Hinshaw et al. (2003b) (5) Hinshaw et al. (2003a) (6) Jarosik et al. (2003b) (7) Kogut 
et al. (2003) (8) Komatsu et al. (2003) (9) Page et al. (2003a) (10) Page et al. (2003b) (11) Peiris et al. (2003) (12) 
Verde et al. (2003). 
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the most difficult to distinguish from true sky signals. The Observatory was designed to minimize all thermal and 
voltage variations and all susceptibilities to these variations, especially at the spin period, as discussed in §2 and by 
Bennett et al. (2003a). In addition, high precision temperature monitors on the Observatory provide the data needed 
to verify that systematic errors from thermal variations are negligible. Jarosik et al. (2003b) report that in-flight spin- 
synchronous effects from the radiometers are < 0.17 fjK rms in the time-ordered-data (TOD), based on flight thermal 
variations multiplied by upper limits on component susceptibilities measured in ground testing. Analysis of flight data 
without use of characterizations derived from ground-based testing give < 0.14 pK rms from all sources (not just the 
radiometers). This is a factor of > 50 times smaller than the requirement that was set in the mission’s systematic error 
budget. Thus, no corrections to the first year WMAPdata are required for spin-synchronous systematic errors . 

The core of the processing pipeline calibrates the data and converts the differential temperatures into maps. The 
data are calibrated based on the Earth-velocity modulation of the CMB dipole. A gain model of the radiometers 
was derived and fit by Jarosik et al. (2003b). The model is based on the constancy of the dipole signal on the sky, 
the measured physical temperature of the front-end radiometer components, and on the time-averaged RF-bias (total 
power) of the radiometer outputs. This relatively simple model closely matches the gains derived from the hourly 
measurements of the amplitude of the dipole and is used in WMAP data processing. Calibration is achieved within 
0.5% accuracy, dominated by the statistical uncertainty in the absolute calibration. 

Low levels of 1 // noise create stripes in the maps that affect the angular power spectrum and other statistics 
derived from the maps. A pre-whitening filter is applied to the TOD to minimize these artifacts. An estimate of the 
magnitude of the striping is given in Hinshaw et al. (2003b) for the maps, and by Hinshaw et al. (2003a) for the power 
spectrum. 

The differential temperature data are formed into maps based on the technique introduced by Wright et al. 
(1996b). HEALPix 11 is used to define map pixels on the sky in Galactic coordinates. Various levels of resolution 
are specified by a “resolution level” with an integer (r = 0, 1,2, ...). With N S jde = 2 r , the number of pixels in the map 
is Np ix = 12 N s 2 ide . The area per pixel is fi p ix = 47 t / A f pix and the separation between pixel centers is 0 pix = fiV 2 . For ex- 
ample, HEALPix resolution level r = 9 (used in WMAP map-making) corresponds with N S jd e = 512, N P j X = 3 145728, 
ftpix = 3.99 x 10- 6 sr, and 0 P j X = 0? 1 15 = 6.87 arc-min. 

WMAP observes the sky convolved with the beam pattern. This is equivalent to the the spatial transform of the 
sky multiplied by the instrument’s “window function.” The beam patterns are measured in-flight from observations of 
Jupiter (Page et al. 2003a). Uncertainties in our knowledge of the beam pattern, although small, are a significant source 
of uncertainty for WMAP since they imply imperfect knowledge of the window function. A small difference between 
the A-side and B-side optical losses was derived based on dipole observations and corrected in the processing. Far 
sidelobes of the beam patterns, determined by ground measurements and in-flight using the Moon, have been carefully 
examined (Barnes et al. 2003). A small far-sidelobe correction is applied only to the K-band map. We now describe 
the maps. 


4. THE MAPS 

We combine the radiometer results within each band and present the five full sky maps at effective CMB 
anisotropy frequencies 23, 33, 41, 61, and 94 GHz in Figures 2a, 2b, 2c, 2d, 2e. The maps are shown in the Mollweide 
projection in units of CMB thermodynamic temperature. The number of independent observations that contribute to 


1 1 http://www.eso.org/science/healpix/ 
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Fig. 2a. — The K-band full sky map in Galactic coordinates, shown in a Mollweide projection. 
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Fig. 2b. — Ka-band 
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Fig. 2c. — Q-band 
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each pixel form the sky pattern in Figure 3. Figure 4 provides an overall guide to some of the more prominent features 
of the maps as well as point sources detected by WMAP, as described by Bennett et al. (2003b). 

Figure 5 shows the K-band and Ka-band maps, with the Ka-band map smoothed to K-band resolution. Note 
both the significant decrease in Galactic signal from K-band to Ka-band and the high Galactic latitude similarities of 
the CMR between the maps. Likewise, Figure 6 shows the Q-band, V-band, and W-band maps with the latter two 
smoothed to Q-band resolution. Higher Galactic contamination in Q-band is apparent. Both Figures 5 and 6 highlight 
the consistency of the high Galactic latitude CMB anisotropy pattern from band to band. 

Comparisons of data between WMAP radiometers, and between WMAP and COBE , are important indicators of 
systematic error levels. Figure 7 illustrates the enormous improvement in angular resolution from COBE to WMAP. 
Features in the maps appear to be generally consistent, but the consistency is better addressed by a more direct compar- 
ison. To do this we take a combination of the WMAP Q-band and V-band maps and smooth it to mimic a COBF-DMR 
53 GHz map (see Figure 8). We then examine the difference between the COBE and WMAP pseudo-map. Figure 9 
shows the difference map along with a simulated map of the noise. With the exception of a feature in the Galactic 
plane, the agreement is clearly at the noise level. The Galactic plane feature is likely to be a result of the spectral index 
uncertainty of combining the Q-band and V-band maps to make a 53 GHz equivalent map. 


5. FOREGROUND ANALYSES 

An understanding of diffuse Galactic emission and extragalactic point sources is necessary for CMB analyses. 
The WMAP mission carries radiometers at five frequencies for the purpose of separating the CMB anisotropy from 
foreground emission based on their different spectra. Figure 10 illustrates the spectral difference between the CMB 
and foregrounds. The WMAP bands were selected to be near the frequency where the ratio of the CMB anisotropy to 
the contaminating foreground is at a maximum. 


5.1. Masks 

For CMB analyses it is necessary to mask out regions of bright foreground emission. Bennett et al. (2003c) 
present a recipe for foreground masks based on K-band temperature levels. Since foreground contamination is most 
severe in K-band, it is used as the best tracer of contamination. The contamination morphology is similar enough 
across all five WMAP bands that masks based on the temperature levels in other bands would be redundant and 
unnecessary. Standard names are given for the mask levels. For example, the KpO mask cuts 21.4% of sky pixels 
while the Kp2 mask cuts 13.1%. See Bennett et al. (2003c) for further detail. An extragalactic point source mask is 
also constructed based on selections from source catalogs. An additional 2% of pixels are masked due to these ~ 700 
sources. 


5.2. Diffuse Galactic emission 

Beyond the use of masks, one technique for reducing the level of foreground contamination is to form a linear 
combination of the multifrequency WMAP data that retains unity response for only the emission component with a 
CMB spectrum. This technique was introduced for COBE by Bennett et al. (1992). With five WMAP bands instead 
of the three on COBE , and with a somewhat more elaborate approach for WMAP ; Bennett et al. (2003c) arrive at the 
internal (WMAPdata only) linear combination map seen in Figure 11 of this paper. The foregrounds are removed to 
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Fig. 3. — The number of independent observations per pixel in Galactic coordinates. The number of observations 
is greatest at the ecliptic poles and in rings around the ecliptic poles with diameters corresponding to the separation 
angle of the two optical boresight directions (approximately 141 °). The observations are the most sparse in the ecliptic 
plane. Small area cuts are apparent where Mars, Saturn, Jupiter, Uranus, and Neptune data are masked so as not to 
contaminate CMB analyses. Jupiter data are used for beam mapping. The histogram of the sky sampling shows the 
departures from uniform sky coverage. 
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Fig. 4. — A guide to the microwave sky for reference. This picture shows the large-scale emission from the Milky Way 
galaxy, including some of its notable components such as the Cygnus complex, the North Polar Spur, the Gum region, 
etc. The small circles show positions of the microwave point sources detected by WMAP (Bennett et al. 2003c). The 
brighter sources are labeled for reference. 
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Fig. 5. — A comparison of the K-band map with the Ka-band map smoothed to K-band resolution, both in thermody- 
namic temperature, shows the dramatic reduction of Galactic contamination with increased frequency. The comparison 
also shows the similarity of the CMB fluctuation features at high Galactic latitude. 
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Fig. 6. — A comparison of the Q-band, V-band, and W-band maps. All three maps are smoothed to Q-band resolu- 
tion and are in thermodynamic temperature. The reduction of Galactic contamination relative to K-band and Ka-band 
(Figure 5) is apparent. The maps show that the constant features across bands are CMB anisotropy while the thermo- 
dynamic temperature of the foregrounds depends on the band (frequency). 
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Fig. 7. — A comparison of the COBE 90 GHz map (Bennett et al. 1996) with the W-band WMAP map. The WMAP 
map has 30 times finer resolution than the COBE map. 
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Fig. 8. — The COjBE-DMR 53 GHz map (Bennett et al. 1996) is shown along with a map made with a linear combi- 
nation of the Q-band and V-band WMAP maps to mimic a 53 GHz map. Note the strong similarity of the maps. 
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Fig. 9. — The difference map is shown between the COBE-DMR 53 GHz map and the combination Q-band/V-band 
maps from Figure 8. This is compared with a map of the noise level. The maps are consistent with one another with 
the exception of a feature in the galatic plane. This discrepancy is likely to be due to a spectral index that is sufficiently 
different from the assumed CMB spectrum used to combine the WMAP Q-band and V-band maps to mimic a 53 GHz 
map. 



Fig. 10. — False color images represent the spectral information from mutliple WMAP bands. Q-band is red, V-band 
is green, and W-band is blue. A CMB thermodynamic spectrum is grey, (top) A three color combination image from 
the Q-, V-, and W-band maps. The dipole and high Galactic latitude anisotropy are seen. ( bottom ) A similar false 
color image but with the dipole subtracted. 


- 21 - 


a remarkable degree; however, the statistics of this internal linear combination map are complex and inappropriate for 
most CMB analyses. 

Below, we use the notation convention that flux density is S ~ and antenna temperature is 7 a ~ where the 
spectral indices are related by /? = a -2. In general, the CMB is expressed in terms of thermodynamic temperature, 
while Galactic and extragalactic foregrounds are expressed in antenna temperature. Thermodynamic temperature 
differences are given by AT = AT A [(e*- l) 2 /jt 2 e*], where x = hvjkT^ h is the Planck constant, v is frequency, k is 
the Boltzmann constant, and T 0 = 2.725 K is the CMB temperature (Mather et al. 1999). Values of AT /A T A for the 
WMAP bands are given by Jarosik et al. (2003b) and can be found in Table 1 . 

Bennett et al. (2003c) identify the amplitudes and spectral indices of the individual emission components. A 
maximum entropy method (MEM) approach is adopted where priors are used for component amplitudes and spectral 
indices, except for free-free emission, which has a fixed spectral index (-2.15 in the WMAP bands). An iterative 
fit is performed, where the pixel-by-pixel amplitudes are updated in accordance with the MEM residuals until low 
(< 1 %) residuals are achieved. The process results in a map of each emission component for each of the five WMAP 
bands. The derived maps of thermal emission from dust give a uniform spectral index across the sky of fid ~ 2.2. The 
derived map of free-free emission is reasonable given the amplitude and morphology of Ha measurements. The other 
radio component fit should include the combined emission of synchrotron and spinning dust. It shows the synchrotron 
spectrum steepening with increasing frequency, as would be expected for a spectral break due to synchrotron losses 
at « 20 GHz. There is no indication of the less steep or flattening spectral index that would result from spinning dust 
emission. Thus spinning dust emission is limited to < 5% of the total Ka-band foreground. Reports of dust-correlated 
microwave emission from COBE data analyses are understood as an admixture of the fraction of synchrotron emission 
(with (3 « -3) that is traced by a dust template, and thermal dust emission (/? « 2.2), giving a combined spectral index 
of /? « -2.2 between the WMAP Ka-band and V-band, approximating the COBE 3 1 GHz and 53 GHz bands. 

While the MEM method is useful for understanding the nature of the foreground emission components, these 
results can not be directly used in CMB analyses due to the complex noise properties that result from the MEM 
process and its simultaneous use of multifrequency maps. For the CMB analyses we use a mask (usually Kp2) to 
exclude pixels where the Galactic emission is strong, combined with template fitting (using external data only) where 
the foregrounds can be adequately corrected. This approach does not complicate the noise properties of the maps. 
Bennett et al. (2003c) describe the template fitting in detail. Thermal dust emission has been mapped over the full 
sky in several infrared bands, most notably by the COBE and IRAS missions. A full sky template is provided by 
Schlegel et al. (1998), and is extrapolated in frequency by Finkbeiner et al. (1999). The mostly synchrotron emission 
map of Haslam et al. (1981) at 408 MHz is used as a radio template. The free- free ionized gas is traced by the Ha 
map assembled by Finkbeiner (2003) from the Wisconsin H-Alpha Mapper (WHAM), the Virginia Tech Spectral- 
Line Survey (VTSS), and the Southern H-Alpha Sky Survey Atlas (SHASSA) (Dennison et al. 1998; Haffner et al. 
2002; Reynolds et al. 2002; Gaustad et al. 2001) . The foreground contamination remaining after subtraction of the 
template-fit component is < 7 /xK rms at Q-band and < 3 fjK rms at both V-band and at W-band for / < 15. 


5.3. Point sources 

A search was made for point sources in the WMAP data. A catalog of 208 detected sources (with 98% reliability) 
is provided by Bennett et al. (2003c). The derived source counts give a power spectrum level of C src = (15 ± 3) x 
10“ 3 /xK 2 sr at Q-band. This is consistent with the level found in the bispectrum analysis of the maps (Komatsu et al. 
2003) and the level found in fits to the map power spectra (Hinshaw et al. 2003b). We have confidence that the point 
source level is understood since it is independently derived using three different methods. 
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Fig. 1 1. — This “internal linear combination” map combines the five band maps in such as way as to maintain unity 
response to the CMB while minimizing foreground foreground contamination. For a more detailed description see 
Bennett et al. (2003c). For the region that covers the full sky outside of the inner Galactic plane, the weights are 
0.109,-0.684,-0.096, 1.92 1,-0.250 for K, Ka, Q, V, and W bands, respectively. Note that there is a chance alignment 
of a particularly warm feature and a cool feature near the Galactic plane. As discussed in Bennett et al. (2003c), the 
noise properties of this map are complex, so it should not generally be used for cosmological analyses. 
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5,4. Sunyaev-Zeldovich (SZ) effect 

Hot gas in clusters of galaxies imparts energy to the CMB photons, causing a temperature decrement in the 
WMAP bands (the Sunyaev-Zeldovich Effect). The Coma cluster is expected to have the most pronounced signature. 
For the highest resolution maps, Bennett et al. (2003c) get -0.34 ±0.18 mK in W-band and -0.24 ± 0. 1 8 mK in V- 
band in the direction of Coma. Use of the XBACS catalog of X-ray clusters as a template results in -0.36 ± 0. 14. This 
verifies that the Sunyaev-Zeldovich Effect is barely detectable in even a matched search of the WMAP sky maps and 
therefore it is not a significant “contaminant” to the WMAPdata. 


6. LIMITS ON NON-GAUSSIANITY 

Maps of the sky are the most complete and compact representation of the CMB anisotropy. Cosmological analyses 
are based on statistical properties of the maps with the power spectrum as one of the most commonly derived statistics. 
The power spectrum is a complete representation of the data only if the CMB anisotropy is Gaussian. Also, the most 
common cosmological models predict that the CMB anisotropy should be consistent with a Gaussian random field (at 
least at levels that are currently possible to measure). Therefore, we test the Gaussianity of the anisotropy, both to 
interpret the power spectrum (and other statistical derivatives of the maps) and to test cosmological paradigms. 

There is no single best test for Gaussianity. Specific tests can be more or less sensitive to different assumed forms 
of non-Gaussianity. Komatsu et al. (2003), in a companion paper, search for non-Gaussianity in the WMAP CMB 
anisotropy maps using Minkowski functionals and a bispectrum estimator. 

Minkowski functionals (Minkowski 1903; Gott et al. 1990) quantify topological aspects of the CMB maps. 
Anisotropy is examined via contours at different temperature levels, and tha number and areas of regions enclosed 
by these contours are computed. Three Minkowski functionals are the area represented by hot and cold spots, the 
contour length around these areas, and the difference between the number of these areas (the “genus”). 

It is widely believed that the CMB anisotropy arises from Gaussian linear fluctuations in the gravitational po- 
tential. Komatsu & Spergel (2001) suggest that non-Gaussian anisotropy be considered as the next, quadratic, term 
of a curvature perturbation Taylor expansion. The coefficient of the quadratic (mean-square) non-Gaussian term is a 
non-linear coupling parameter, fy L . Since the CMB bispectrum measures the phase correlations of the anisotropy, it 
can be used to solve for fyi. The Minkowski functional results can also be expressed in terms of fyi. 

For the Minkowski functionals, Komatsu et al. (2003) find < 139 (95% CL). From the bispectrum, Komatsu 
et al. (2003) find -5 8 < fyL < 134 (95% CL). The two results are consistent. The CMB anisotropy is thus demonstrated 
to follow Gaussian statistics. This is a significant result for models of the early universe. It also means that we can 
construct and interpret CMB statistics (e.g. the angular power spectrum) from the maps in a straightforward manner. 


7. MULTIPOLES 

The temperature anisotropy, T^n), is naturally expanded in a spherical harmonic basis, F/ m , as 

T( n ) = Y^a,Mn). ( 1 ) 

The angular power spectrum, C/, is a cosmological ensemble average given by 

c, = (M 2 > 


( 2 ) 
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and is observed for our actual sky as 

0) 

m 

Assuming random phases, the temperature anisotropy for each multipole moment, A 7}, can be associated with the 
angular spectrum, C/, as 

A7] = yj Cf c> ' 1(1 + 1)/2tt. (4) 

The correlation function is 

C(6) = 


where P/ is the Legendre polynomial or order l. 

In practice, an instrument adds noise and spatially filters the sky signal due to the beam pattern and any other ex- 
perimental limitations on the sampling of all angular scales. The experimental transfer function, bu for each multipole 
moment depends upon the specific beam pattern of the experiment. WMAP samples all angular scales limited only by 
the beam pattern. The window function for the signal power between channels i and j is 

w) J =b i ! b J ,p}, (7) 

where pj is the pixel window function. The angular power spectrum then becomes 

C^d/bibj+NiSij ( 8 ) 

where N\ is the power spectrum noise that results from the instrument noise, which is assumed to be uncorrelated 
between channels i and j. 

The use of a sky mask for foreground suppression breaks the orthogonality of the spherical harmonics on the sky 
and leads to mode coupling. Hinshaw et al. (2003b) discuss how WMAP handles this, and other complexities. 


(5) 

1 m=-l 

L£(2/+1)C,PK«w0) (6) 


7.1. 1 = 1 dipole 

COBE determined the dipole amplitude is 3.353 ± 0.024 mK in the direction (/, b) = (264.26° ± 0.33°, 48.22° ± 
0.13°), where / is Galactic longitude and b is Galactic latitude (Bennett et al. 1996). This dipole was subtracted from 
the WMAP data during processing. Examination of the WMAP maps allow for the determination of a residual dipole, 
and thus an improvement over the COBE value. Note that this does not have any effect on WMAP calibration, which 
is based on the Earth’s velocity modulation of the dipole, and not on the dipole itself. The WMAP -determined dipole 
is 3.346 ± 0.017 mK in the direction (/, b)= (263P85 ± 0? 1, 48?25 ± 0?04). The uncertainty of the dipole amplitude is 
limited by the WMAP 0.5% calibration uncertainty, which will improve with time. 


7.2. I = 2 quadrupole 

The quadrupole is the l ~ 2 term of the spectrum A 7} 2 = /(/ + 1)C//27 t, i.e. A7]^ 2 - (3/7r)C/ = 2. Alternately, the 
quadrupole amplitude can be expressed as Q rms = >/ (5 /47r) C /=2 = i /5jTZ A7^ 2 - The 4-year COBE quadrupole is 



-25- 


Qrms = IOI4 /xK with the peak of the likelihood in the range 6.9 /xK < Q rms < 10 jiK, as shown in Figure 1 of Hinshaw 
etal. (1996a). 

The WMAP quadrupole, Q rms = 8 ± 2 /xK or AT 2 2 = 1 54 ± 70 /xK 2 , is consistent with COBE but with tighter limits 
because of better measurements and understanding of foregrounds. The foreground level is still the leading uncertainty. 
(The small kinematic quadrupole is not removed from the maps nor accounted for in this analysis.) The quadrupole 
value is low compared with values predicted by ACDM models that fit the rest of the power spectrum. ACDM models, 
in particular, tend to predict relatively high quadrupole values due to the enhanced, A-driven, integrated Sachs- Wolfe 
effect. 


7.3. n-poles 

A central part of the task of computing multipole information from the maps is the evaluation and propagation 
of errors and uncertainties. This largely involves arriving at an adequate representation of the Fisher matrix, which is 
the inverse covariance matrix of the data. The Fisher matrix must take into account mode-coupling from the sky cut, 
beam (window function) uncertainties, and noise properties. The fact that the WMAP data use a nearly azimuthally 
symmetric cut, and have a nearly diagonal pixel-pixel covariance, greatly simplifies the evaluation of the Fisher matrix. 

Two approaches to computing the angular power spectrum have been used by Hinshaw et al. (2003b): a quadratic 
estimation based on Hivon et al. (2002); and a maximum likelihood estimate based on Oh et al. (1999). The quadratic 
estimator is used in the final WMAP spectrum analyses, while the maximum likelihood technique is used as a cross- 
check. 

The K-band and Ka-band beam sizes are large enough that these bands are not used for CMB analysis since 
they have the most foreground contamination and probe the region in /-space that is cosmic variance limited by the 
measurements at the other bands. These bands are invaluable, however, as monitors of Galactic emission. The two 
Q-band, two V-band, and four W-band differencing assemblies are the source of the prime CMB data. The matrix of 
auto- and cross-correlations between the eight Q-, V-, and W-band differencing assemblies has eight diagonal (auto- 
correlations) and 28 unique off-diagonal elements (cross-correlations). Since auto-correlations are difficult to assess 
due to the noise bias (see equation 8), we included only the 28 unique off-diagonal (cross-correlations) in the WMAP 
power spectrum analysis. In dropping the auto-correlations, each of which has twice the noise variance of a cross- 
correlation, we lose only 1 - i/56/(56+8) = 6% of the ideally achievable signal-to-noise ratio. In Hinshaw et al. 
(2003b), we show that the power spectrum computed from the auto-correlation data is consistent with the angular 
power spectrum from the cross-correlation data. We anticipate using the auto-correlation data in future analyses. 

The cross-power spectrum from the 28 pairs is considered in four 1 ranges. For l < 100 we use uniform pixel 
weighting of only V- and W-band data. This reduces the Galactic contamination where measurement errors are well 
below the cosmic variance. For 100 < l < 200 we use uniform pixel weighting of the combined 28 pairs. For 200 < 
/ < 450 all 28 cross-power pairs are used with a transitional pixel weighting. For / > 450 all 28 pairs are used with 
inverse noise weighting. Our Monte Carlo simulations show that this is a nearly optimal scheme. 

The angular power spectrum is shown for the WMAP data in Figure 12. The WMAP power spectrum agrees 
closely with COBE at the largest angular scales, and with CBI and ACBAR at the finer angular scales. We highlight the 
CBI and ACBAR results because they are a useful complement to WMAP at the smaller angular scales. The acoustic 
pattern is obvious. Page et al. (2003b) find that the first acoustic peak is A 7) = 74.7 ± 0.5 jjK at l = 220.1 ± 0.8. The 
trough following this peak is 41 .0 ± 0.5 /xK at / = 41 1.7 ± 3.5 and the second peak is 48.8 ± 0.9 /xK at / = 546 ± 10. 

ACDM models predict enhanced large angle power due to the integrated Sachs-Wolfe effect. The WMAP and 
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Fig. 12. — The WMAP angular power spectrum, (top:) The WMAP temperature (TT) results are consistent with the 
ACBAR and CBI measurements, as shown. The TT angular power spectrum is now highly constrained. Our best fit 
running index ACDM model is shown. The grey band represents the cosmic variance expected for that model. The 
quadrupole has a surprisingly low amplitude. Also, there are excursions from a smooth spectrum (e.g., at t « 40 and 
t k, 210) that are only slightly larger than expected statistically. While intriguing, they may result from a combination 
of cosmic variance, subdominant astrophysical processes, and small effects from approximations made for this first 
year data analysis (Hinshaw et al. 2003b). We do not attach cosmological significance to them at present. More 
integration time and more detailed analyses are needed, (bottom:) The temperature-polarization (TE) cross-power 
spectrum, (/ + l)Q/2n. (Note that this is not multiplied by the additional factor of /.) The peak in the TE spectrum 
near l ~ 300 is out of phase with the TT power spectrum, as predicted for adiabatic initial conditions. The antipeak in 
the TE spectrum near / ^ 150 is evidence for superhorizon modes at decoupling, as predicted by inflationary models. 
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COBE data, on the other hand, have the opposite trend. The conflict is also seen clearly in the correlation function, 
C(0 ), shown in Figure 13. The WMAP correlation function is computed on the KpO cut sky, and the COBE correlation 
function is computed on the “custom cut” sky (Bennett et al. 1996). The best-fit ACDM model is shown with a grey 
band indicating one standard deviation as determined by Monte Carlo simulations. 

The lower two plots in Figure 13 display the correlation function of the difference between the COBE-DMR 
and WMAP maps with a \b\ = 10° Galactic plane cut. A synthesis of the WMAP Q- and V-band maps was made 
to approximate a 53 GHz map to compare with the CORE-DMR 53 GHz map. The COBE-DMR 90 GHz map is 
compared directly, without corrections, to the WMAP W-band map. These plots emphasize the consistency of the 
WMAP and COBE measurements. The slightly higher than expected deviations at 53 GHz are likely to be due to 
Galactic contamination, arising from outside the cut regions and from the construction of the synthesized WMAP 53 
GHz map. 

Since a correlation plot emphasizes the very low / power, the discrepancy between the model and WMAP data 
illustrates that there is surprisingly little anisotropy power in the WMAP (and COBE) maps at large separation. This 
result is generically true for ACDM models, independent of the exact parameters of the model. There is very little 
large scale CMB anisotropy power in our sky . This fact, first seen by COBE is confirmed by WMAP. The probability 
of so little C($ > 60°) anisotropy power is ~ 2 x 10" 3 , given the best-fit ACDM model. The lack of power is seen both 
in C(0) and the behavior of the low order (/ = 2, 3, 4, and 5) multipoles. 


8. CMB POLARIZATION & THE DETECTION OF REIONIZATION 

Each differencing assembly measures the sky in two orthogonal linear polarizations. As the Observatory spins, 
precesses, and orbits the Sun, the instrument observes the sky over a range of polarization angles. The range of angles 
observed is neither complete nor uniform, but it is sufficient to provide valuable new CMB polarization results. 

We express polarization in term of the standard /, Q, and U Stokes parameters, in units of thermodynamic 
temperature. Polarization coordinates are not rotationally invariant so a coordinate system is defined by Kogut et al. 
(20Q3) for expressing information. By appropriate summing and differencing of the time-ordered-data from the pairs 
of corresponding radiometers, we iteratively make maps in unpolarized intensity, /, and in each polarized component, 
Q and U (Hinshaw et al. 2003b). The Q and U maps not only have much lower signal-to-noise than the / maps, but 
they are more susceptible to systematic errors since the differencing between radiometer outputs occurs on the ground, 
not in the front-end of the radiometers. 

The U polarization maps are expected to be most affected by systematic errors. This is due to the orientation of 
the observed polarization angles relative to the spin axis of the Observatory, causing systematics not fixed on the sky 
to preferentially go into the U map. 

Taurus A, an extremely strong source in the WMAP bands, is polarized and its polarization amplitude and di- 
rection are known (Flett & Henderson 1979). Polarization measurements from WMAP are consistent with previous 
observations of Taurus A, providing a useful check on the operation of the WMAP hardware and software. However, 
the polarization systematic errors have not yet been fully quantified. We chose to release the / data and results as soon 
as possible, rather than postpone their publication and release until the systematic measurement errors of the Q and 
U data are fully assessed. All data, including the Q and U maps will be released when the characterization of their 
instrumental signatures are complete. Results of temperature-polarization correlations between the / and Q maps, i.e. 
TE correlations, are much less sensitive to systematic effects than polarization signals alone, so we are able to report 
those results. 
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Fig. 13. — (top): CMB temperature correlation function of the WMAP and COBE data. The WMAP correlation 
function is computed using the KpO cut sky, and the COBE correlation function is computed using the “custom cut” 
sky. The running index ACDM model that is fit to the power spectrum is shown with a Monte Carlo determined grey 
band indicating one standard deviation. The model is, overall, an excellent fit to the WMAPpower spectrum. However, 
a correlation plot emphasizes the low / power. The discrepancy between the model and data illustrates that there is 
surprisingly little anisotropy power in the WMAP and COBE maps at large angles, (bottom): The lower two plots 
display the correlation function of the difference between the COBB-DMR and WMAP maps with a \b\ = 10° Galactic 
plane cut. A synthesis of the WMAPQ and V band maps was made to approximate a 53 GHz-like map to compare 
with the CQBE-DMR 53 GHz map. The COBE-DMR 90 GHz map is compared directly, without corrections, to the 
WMAP W-band map. These plot emphasize the consistency of the WMAP and COBE measurements. The slightly 
higher than expected deviations at 53 GHz are likely to be due to Galactic contamination, arising from outside the cut 
regions and from the construction of the synthesized WMAP 53 GHz map. 
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Figure 12 shows the temperature-polarization (TE) cross power spectrum. Kogut et al. (2003) report a detection of 
TE correlations on both large angular scales (from reionization) and on small scales (from the adiabatic fluctuations). 
The TE power spectrum, shown in Figure 12, is discussed in detail by Kogut et al. (2003). In the TE angular power 
spectrum the antipeak is -35 ±9 fiK 2 at /= 137 ±9 and the peak is 105 ± 18 fjK 2 at l = 329 ± 19 (Page et al. 2003b). 

The detection of the reionization of the universe corresponds to an optical depth r = 0. 17 ± 0.04 (0.09 < r < 0.28 
at 95% confidence). Although WMAP measures the integrated optical depth, the epoch and redshift of reionization 
can be derived firom the integral optical depth within the context of a model of the reionization process. A single 
instantaneous step function from r = 0 to a steady fixed value of r is physically unlikely. For more likely models, Kogut 
et al. (2003) conclude that the redshift of reionization is z T = 20^° (95% CL), corresponding to an epoch of reionization 
of t r = 180^|q° Myr (95% CL) after the Big Bang. Cen (2002) presents a detailed model of a reionization process and 
predicts an integral value of the optical depth, r = 0. 10 ± 0.03. This is low compared with the newly measured value, 
but would be higher and consistent if the assumed cosmological parameters in the model were adjusted to the values 
of the new WMAP best-fit parameters (see §9.1). 

The measured optical depth means that reionization suppressed the acoustic peak amplitudes by l - e~ 2r « 30%. 
While accounted for in our model fits (§9, below), this suppression was not accounted for in previous CMB parameter 
determinations. The anticorrelations observed in the TE power spectrum directly imply superhorizon fluctuations, a 
new result in support of inflation-like theories, as discussed in §9.3. 


9. COSMOLOGICAL INTERPRETATION 

In this section we summarize the cosmological intepretation of WMAP first year results, which are discussed in 
more detail by Spergel et al. (2003), Peiris et al. (2003), Page et al. (2003b), and Kogut et al. (2003). The methodology 
used in the model fits is described by Verde et al. (2003). 

Spergel et al. (2003) show that a cosmological model with a flat universe, seeded with a scale-invariant spectrum 
of adiabatic Gaussian fluctuations, with reionization, is an acceptable fit not only to the WMAP data but also to a host 
of astronomical data. These data are: smaller angular scale CMB anisotropy data from ACBAR (Kuo et al. 2002) and 
CBI (Pearson et al. 2002); the HST key project value of Ho (Freedman et al. 2001); the accelerating Universe seen in 
Type la SNe (Riess et al. 2001); the shape and amplitude of the large scale structure seen in clusters and superclusters 
of galaxies (Percival et al. 2001 ; Verde et al. 2003); and the linear matter power spectrum seen in the Lyman a forest 
(Croft et al. 2002). There has been mounting evidence in the direction of this model for years (Peebles 1984; Bahcall 
et al. 1999). The optical depth since reionization is a new, but not surprising component of the model. The WMAP 
data establish this model as the standard model of cosmology by testing the key assumptions of the model and by 
enabling a precise determination of its parameters. 

The WMAP data test several of the key tenets of the standard model. The WMAP detection of temperature- 
polarization correlations (Kogut et al. 2003) and the clear detection of acoustic peaks (Page et al. 2003b) implies that 
the primordial fluctuations were primarily adiabatic: the primordial ratio of dark matter/photons and the primordial 
ratio of baryons/photons do not vary spatially. The analysis of the WMAP temperature data demonstrates Gaussianity 
(Komatsu et al. 2003). The WMAP data, when combined with any one of the following three external data sets: HST 
Key Project measurement of H 0 (Freedman et al. 200 1 ), the 2dFGRS measurement of the matter density (Percival et al. 
2001; Verde et al. 2003) or the Type la supernova measurements (Riess et al. 2001) imply that the radius of curvature 
of the universe,/? = c//o 1 |l-n mf |“ 1/2 , must be very large, fl to t= 1.02±0.02 and 0.99 < (l tot < 1.05 (95% CL). These 
measurements also require that the dark energy be the dominant constituent of fl tot . The WMAP data alone rule out 
the standard = 1 CDM model by > la . 
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9,1. Best Fit Cosmological Model 

While an acceptable fit, the model described above is not our best fit model. In the discussion below we con- 
centrate on our best fit model, which adds a scale-dependent primordial spectral index. This cosmological model is a 
flat universe with a baryon fraction of ft* = 0.044 ± 0.004, a matter fraction of ft m = 0.27 ± 0.04, and a dark energy 
fraction of ft a = 0.73 ± 0.04, seeded with a scale-dependent spectrum of adiabatic Gaussian fluctuations. This model 
has C/= io = 46.0 /iK 2 , consistent with the COBE measurement of Q = w = 44.4 /xK 2 . 

The WMAP data alone enable accurate determinations of many of the key cosmological parameters (Spergel 
et al. 2003). But a combination of the WMAP data with the COBE determination of the CMB temperature (Mather 
et al. 1999), the CBI (Pearson et al. 2002) and the ACBAR (Kuo et al. 2002) CMB measurements, and the 2dFGRS 
survey determination of the power spectrum of the local galaxy fluctuations (Percival et al. 2001), yields the best 
fit cosmological parameters listed in Table 3. Verde et al. (2003) describes our methodology for determining these 
parameters and Spergel et al. (2003) describes the best fit models for different combinations of data sets. 


A power spectrum of primordial mass fluctuations with a scale invariant spectral index is given by P(k) = Ak nx 
with n s - d In P/ dink. Inflationary models predict a running spectral index (Kosowsky & Turner 1995), and our best 
fit model uses a power spectrum of primordial mass fluctuations with a scale-dependent spectral index: 


P(k)^P(k o) 



n s (k 0 y+( \/l){dn s jd In k) ln(fc/*o) 


As in the scale-independent case, we define 


n s (k) = 


dlnP 
dink ’ 


so 


n s (k) = n s (ko) + 


dn s 

dink 


In 



(9) 

GO) 
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(with d 2 n s /d\nk 2 = 0). The definition for n s used here includes a factor of (1/2) difference from the Kosowsky & 
Turner (1995) definition. Peiris et al. (2003) explore the implications of this running spectral index for inflation. 
The best fit values of A, n s and dn s /d\nk are in Table 3 for ko = 0.05 Mpc -1 . A is the normalization parameter in 
CMBFAST version 4. 1 with option UNNORM. The amplitude of curvature fluctuations at the horizon crossing is 
| A/?(£o)| 2 = 2.95 x 10~ 9 A. We discuss the implications of the measured values of these parameters in §9.3 and in Peiris 
et al. (2003). 

The WMAP data constrains the properties of both the dark matter and the dark energy in the following ways: 


(a) The WMAP detection of reionization at z ~ 20 is incompatible with the presence of significant warm dark matter 
density. Since the warm dark matter moves too fast to cluster in small objects, the first objects do not form in 
this scenario until z ~ 8 (Barkana et al. 2001). 

(b) The running spectral index implies a lower amplitude for mass fluctuations on the dwarf galaxy scale. Dark 
matter simulations of models (Ricotti 2002) find that the dark matter mass profiles depend upon the spectral 
index on the relevant mass scale. Thus, the shallower spectral index implied by our best fit model may solve the 
CDM dark matter halo profile problem (Moore et al. 1998; Spergel & Steinhardt 2000). 

(c) While the WMAP data alone are compatible with a wide range of possible properties for the dark energy, the 
combination of the WMAP data with either the HST key project measurement of Hq , the 2dFGRS measurements 
of the galaxy power spectrum or the Type la supernova measurements requires that the dark energy be 73% of 
the total density of the Universe, and that the equation of state of the dark energy satisfy w < -0.78 (95% CL). 
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9.2. Best Fit Parameters 

WMAP s measurements of the baryon density, Hubble constant, and age of the universe strengthen the cosmic 
consistency that underlies the Big Bang model. 

ATOMIC DENSITY (Vtbh 1 ): WMAP measures the atomic density at recombination to an accuracy of 4% through 
the shape of the angular power spectrum, and particularly through the ratio of the heights of the first to second peak 
(Page et al. 2003b; Spergel et al. 2003). Our best fit value is Q> b h 2 = 0.0224 ± 0.0009. The baryon density is also probed 
via abundance measurements of [D]/[H] (O’Meara et al. 2001; Pettini & Bowen 2001; D’Odorico et al. 2001). It is 
impressive that £l b h 2 is the same at z = 1 089 as measured via the CMB as it is at z = 1 0 9 from Big Bang nucleosynthesis. 
Thus we find cosmic consistency of the baryon density throughout cosmic time and measurement technique. 

HUBBLE CONSTANT (Ho): The WMAP measurements of the age and Pt m h 2 yield a measurement of Ho = 
71±3 km s -1 Mpc -1 that is remarkably consistent with the HST Key Project value of Ho = 72 ± 3 ± 7 km s _1 Mpc -1 
(Freedman et al. 2001), but with smaller uncertainty. Recent measurements of the Hubble constant from gravitational 
lens timing and the Sunyaev-Zeldovich effect yield independent estimates that are generally consistent, but with larger 
uncertainties at present. Through a variety of measurement techniques that sample different cosmic times and distances 
we find cosmic consistency on // 0 . 

AGE OF THE UNIVERSE (to): The first acoustic peak in the CMB power spectrum represents a known acoustic 
size ( r s = 147 ± 2 Mpc) at a known redshift (Zdec = 1089 ± 1). From these, WMAP measures the age of the universe 
(r 0 = 13.7 ± 0.2 Gyr) to an accuracy of ^ 1% by determining the CMB light travel time over the distance determined 
by the decoupling surface (d A = 14.0^ Gpc) and the geometry of the universe (i.e., flat). The age of the universe is 
also estimated via stars in three ways: 

(i) the main sequence turn-off in globular clusters yielding a cluster age of 12 ± 1 Gyr (Reid 1997); 

(ii) the temperature of the coldest white dwarfs in globular clusters yielding a cluster age of 12.7 ± 0.7 Gyr (Hansen 
et al. 2002), and 

(iii) nucleosynthesis age dating yielding an age of 15.6 ± 4.6 Gyr (Cowan et al. 1999). 

These stellar ages are all consistent with age of the universe found by WMAP. 

MATTER DENSITY (Tt m h 2 ): The matter density affects the height and shape of the acoustic peaks. The baryon- 
to-matter ratio determines the amplitude of the acoustic wave signal and the matter-to-radiation ratio determines the 
epoch, z eq , when the energy density of matter equals the energy density of radiation. The amplitude of the early 
integrated Sachs- Wolfe effect signal is sensitive to the matter-radiation equality epoch. From these effects WMAP 
measures the matter density, fl m h 2 ,to an accuracy of ~ 5%. Large scale structure observations measure Q m h through 
the shape of the power spectrum. When combined with estimates of h , this yields Ct m h 2 . Large scale velocity field 
measurements yield , where b is the bias in how the galaxy power spectrum traces the matter power spectrum 

(/gal = b 2 P(k)). Galaxy bispectrum measurements yield b , allowing for estimates of From the galaxy data, Verde 
et al. (2002) find = 0.27 ± 0.06, which is consistent with the WMAP result of = 0.27 ±0.04. 

Cluster lensing observations yield measurements of the total mass in the cluster. X-ray measurements give both 
the baryonic mass and the total mass through the gravitational potential. Sunyaev-Zeldovich effect observations give 
a different determination of the baryonic mass in clusters. The combined X-ray and SZ measurements give a value of 
QkQ” 1 = 0.081^^/r 1 (Grego et al. 2001 ), which is reasonably consistent with 0.17 ± 0.01 from WMAP. 
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9.3. Implications for Inflation 

WMAP data tests several of the key predictions of the inflationary scenario (see Peiris et al. (2003) for further 
discussion): 

(a) Inflation predicts that the universe is flat. As noted in §9.1 and discussed in detail in Spergel et al. (2003), the 
combination of WMAPdata with either // 0 , Type la SNe, or large scale structure data constrains |1 -PL tot \ <0.03. 

(b) Inflation predicts Gaussian random phase fluctuations. Komatsu et al. (2003) shows that the CMB fluctuations 
have no detectable skewness and place strong constraints on primordial non-Gaussianity. Komatsu et al. (2003) 
also shows that the Minkowski functionals of the WMAP data are consistent with the predictions of a model 
with Gaussian random phase fluctuations. 

(c) Inflation predicts fluctuations on scales that appear to be superhorizon scales in a Friedman-Robertson- Walker 
(FRW) cosmology. The WMAP detection of an anti-correlation between polarization and temperature fluc- 
tuations on scales of ~ l°-2° (Kogut et al. 2003) confirms this prediction and rules out subhorizon causal 
mechanisms for generating CMB fluctuations (Peiris et al. 2003). 

(d) Inflation predicts a nearly scale invariant spectrum of fluctuations, as seen by WMAP. 

The WMAP data, in combination with complementary cosmological data, not only test the basic ideas of the 
inflationary scenario but also rule out broad classes of inflationary models, and therefore the data guide us towards a 
specific workable inflationary scenario. The WMAP data place significant constraints on r, the tensor-to-scalar ratio, 
n s , the slope of the scalar fluctuations and dn s /d Ink, the scale dependence of these fluctuations. The addition of an 
admixture of isocurvature modes does not improve the WMAP model fits. 

The best fit model to the combination of the WMAP, ACBAR, CBI, 2dFGRS and the Lyman-a forest data has 
a spectral index that runs from n > 1 on the large scales probed by WMAP to n < 1 on the small scales probed by 
the 2dFGRS and the Lyman-a forest data. Only a handful of inflationary models predict this behavior. The Linde & 
Riotto (1997) hybrid inflationary model is one example. The data, however, do not yet require n > 1 on large scales: 
our best fit model has n s = 1 .03 ± 0.04 at k = 0.002 Mpc~* . 

Our analysis of inflationary models (Peiris et al. 2003) marks the beginning of precision experimental tests of spe- 
cific inflationary models. With the addition of on-going WMAP observations and future improved analyses, WMAP 
will be able to more accurately constrain r and hence n 5 on large scales. When other CMB experiments are calibrated 
directly to the WMAP sky maps, they will provide improved measurements of the temperature angular power spec- 
trum for / > 700. The upcoming release of the Sloan Digital Sky Survey (SDSS) power spectrum will provide an 
improved measurement of the galaxy power spectrum. The SDSS Lyman-a forest data are expected to be a significant 
improvement over the data used in our current analysis. Looking further towards the future, ESA’s PLANCK mission 
will provide improved measurements on the CMB angular power spectrum on smaller angular scales and should be 
able to improve constraints on r. 


10. DATAPRODUCTS 

All of the WMAP data will be released. In addition, several ancillary and analyzed data sets are released. These 
include beam patterns, angular spectra, etc. Some software tools are also provided. An Explanatory Supplement 
provides detailed information about the WMAP in-flight operations and data products (Limon et al. 2003). All WMAP 
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data products are distributed through the Legacy Archive for Microwave Background Data Analysis (LAMBDA) at 
http://iambda.gsfc.nasa.gov. This is a new NASA data center dedicated to the rapidly growing field of microwave 
background data archiving and analysis. 


11. SUMMARY AND CONCLUSIONS 

(1) WMAP has produced high quality full sky maps in five widely separated frequency bands. These maps can 
be used to test cosmological models and serve as the primary legacy of the mission. 

(2) We have characterized and placed stringent limits on systematic measurement errors. The calibration is based 
on the modulation of the CMB dipole, and is accurate to better than 0.5%. 

(3) We have demonstrated the ability to separate the CMB anisotropy from Galactic and extragalactic foregrounds. 
We provide masks for this purpose. In addition, we have produced CMB maps where the Galactic signal is minimized. 

(4) We have a new determination of the dipole. It is 3.346 ± 0.017 mK in the direction (/, b)- (263.°85 ± 0? 1, 
48?25db0?04). 

(5) We have a new determination of the quadrupole amplitude. It is Q rms = 8 ± 2 ixK or A7 2 2 = 154 ±70 /*K 2 . 

(6) We have placed tight new limits on non-Gaussianity of the CMB anisotropy. The coupling coefficient of a 
quadratic non-Gaussian term is limited to -58 < fyi < 134 (95% CL) (Komatsu et al. 2003). 

(7) We have produced an angular power spectrum of the anisotropy with unprecedented accuracy and precision. 
The power spectrum is cosmic variance limited for / < 354 with a signal-to-noise ratio >1 per mode to / = 658. 

(8) We have, for the first time, observed the angular power spectrum of TE temperature-polarization correlations 
with sufficient accuracy and precision to place meaningful limits on cosmology. 

(9) We have detected the epoch of reionization with an optical depth of r = 0. 17 ±0.04. This implies a reionization 
epoch of t r = 180^|o° Myr (95% CL) after the Big Bang at a redshift of z r = 20^° (95% CL) for a range of ionization 
scenarios. This early reionization is incompatible with the presence of a significant warm dark matter density. 

(10) We have fit cosmological parameters to the data. We find results that are consistent with the Big Bang theory 
and inflation. We find that the addition of a running spectral index, while not required, improves the fit at the ~ 2 a 
level. We provide values and uncertainties for a host. of parameters based on this non-power-law inflationary model. 
Our “best” values for cosmic parameters are given in Table 3. 

(11) WMAP continues to collect data and is currently approved for 4 years of operations at L 2 . The additional 
data, and more elaborate analyses, will help to further constrain models. The addition of other continuously im- 
proving CMB and large scale structure observations is essential for progress towards the ultimate goal of a complete 
understanding of the global properties of the universe. 


The WMAP mission is made possible by the support of the Office of Space Sciences at NASA Headquarters and 
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Table 2. Data Flagging Summary 


Category 

K-Band 

Ka-Band 

Q-Band 

V-Band 

W-Band 

Rejected or Lost Data 

Lost Data (%) 

0.27 

0.27 

0.27 

0.27 

0.27 

Spacecraft thermal change (%) 

0.87 

0.87 

0.87 

0.87 

0.87 

Gain or baseline step (%) 

0.00 

0.13 

0.12 

0.00 

0.22 

Total Lost or Bad Data (%) 

1.04 

1.27 

1.26 

1.14 

1.36 

Data Not Used in Maps 

Planet flag (%) 

0.11 

0.11 

0.11 

0.11 

0.11 
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Table 3. “Best” Cosmological Parameters 


Description 

Symbol 

Value 

+ uncertainty 

- uncertainty 

Total density 

Qtot 

1.02 

0.02 

0.02 

Equation of state of quintessence 

W 

< -0.78 

95% CL 

— 

Dark energy density 

Ha 

0.73 

0.04 

0.04 

Baryon density 

n** 2 

0.0224 

0.0009 

0.0009 

Baryon density 

Vlb 

0.044 

0.004 

0.004 

Baryon density (cm" 3 ) 

rib 

2.5 x 10' 7 

0.1 x 10" 7 

0.1 x lO -7 

Matter density 

n m h i 2 

0.135 

0.008 

0.009 

Matter density 

Cl m 

0.27 

0.04 

0.04 

Light neutrino density 

Q u h 2 

< 0.0076 

95% CL 

— 

CMB temperature (K) a 

Tomb 

2.725 

0.002 

0.002 

CMB photon density (cm" 3 ) b 


410.4 

0.9 

0.9 

Baryon-to-photon ratio 

r\ 

6.1 x 10-'° 

0.3 x lO' 10 

0.2 x 10" 10 

B ary on-to -matter ratio 


0.17 

0.01 

0.01 

Fluctuation amplitude in 8/T 1 Mpc spheres 


0.84 

0.04 

0.04 

Low-z cluster abundance scaling 

<* 8^m 5 

0.44 

0.04 

0.05 

Power spectrum normalization (at ko = 0.05 Mpc" 1 

) C ^ 

0.833 

0.086 

0.083 

Scalar spectral index (at ko = 0.05 Mpc _1 ) c 

n s 

0.93 

0.03 

0.03 

Running index slope (at = 0.05 Mpc" 1 ) 0 

dn s /d]nk 

-0.031 

0.016 

0.018 

Tensor-to-scalar ratio (at ko = 0.002 Mpc" 1 ) 

r 

<0.71 

95% CL 

— 

Redshift of decoupling 

Zdec 

1089 

1 

1 

Thickness of decoupling (FWHM) 

A.Zdec 

195 

2 

2 

Hubble constant 

h 

0.71 

0.04 

0.03 

Age of universe (Gyr) 

to 

13.7 

0.2 

0.2 

Age at decoupling (kyr) 

tdec 

379 

8 

7 

Age at reionization (Myr, 95% CL)) 

tr 

180 

220 

80 

Decoupling time interval (kyr) 

Atdec 

118 

3 

2 

Redshift of matter-energy equality 

Zeq 

3233 

194 

210 

Reionization optical depth 

T 

0.17 

0.04 

0.04 

Redshift of reionization (95% CL) 

Zr 

20 

10 

9 

Sound horizon at decoupling (°) 

Oa 

0.598 

0.002 

0.002 

Angular diameter distance to decoupling (Gpc) 

&a 

14.0 

0.2 

0.3 

Acoustic scale d 

Ia 

301 

1 

1 

Sound horizon at decoupling (Mpc) d 

r s 

147 

2 

2 


a from COBE (Mather et al. 1999) 
b derived from COBE (Mather et al. 1999) 
c l e f f « 700 

% = ttO^ 1 6 a = r s d ~ 1 


